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Abstract 

The subject of this paper is the empirically determined anomalous 
secular increases of the astronomical unit, of the order of some cm yr~^, 
and of the eccentricity of the lunar orbit, of the order of 10^^^ yr~^. 
The aim is to find an empirical explanation of both the anomalies, as far 
as their orders of magnitude are concerned. The methods employed are 
working out perturbatively with the Gauss equations the secular effects 
on the semi- major axis a and the eccentricity e of a test particle orbiting 
a central body acted upon by a small anomalous radial acceleration 
A proportional to the radial velocity Vr of the particle-body relative 
motion. The results show that non-vanishing secular variations (d) 
and (e) occur. If the magnitude of the coefficient of proportionality of 
the extra-acceleration is of the same order of magnitude of the Hubble 
parameter Hq = 7.47 x 10~^^ yr~^ at the present epoch, they are 
able to explain both the astrometric anomalies without contradicting 
other existing observational determinations for the Moon and the other 
planets of the solar system. Finally, it is concluded that the extra- 
acceleration might be of cosmological origin, provided that the relative 
radial particle-body motion is accounted for in addition to that due to 
the cosmological expansion only. Anyway, further data analyses should 
confirm or disproof the existence of both the astrometric anomalies as 
genuine physical phenomena. 

Keywords: celestial mechanics; ephemerides; gravitation; Moon; planets and 
satellites: general 
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1 Introduction 



Recently, the main features of the anomalous secular increases of both the 
astronomical unit and the eccentricity e of the lunar orbit have been reviewed 
[1]. While the first effect, obtained by several independent researchers [21 
El m O [1], should be of the order of a few cm yr~^, the second one [6l [7| 
amounts to e = (9 ± 3) x 10~^^ yr"^, according to the latest data analysis 

m- 

Such phenomena attracted the attention of various scientists dealing with 
them in different contexts PIIHlIIIliniinilllllSlIISlIIZlIISlIiaiSniEI]. 
Thus, several more or less sound attempts to find, or to rule out, possible 
explanations [llllSlllSllMlllSllMlllZllMISlEQliniE^ 
EZlEHlESllinKIlllIlllSlllS] for both the anomahes were proposed so far, 
both in terms of standard known gravitational physical phenomena and of 
long-range modified models of gravity. 

Here we propose an empirical formula which is able to accommodate both 
the anomalies, at least as far as their orders of magnitude are concerned. 

2 An anomalous acceleration proportional to the 
radial velocity of the test particle 

Let us assume that, in addition to the usual Newtonian inverse-square law 
for the gravitational acceleration imparted to a test particle by a central 
body orbited by it, there is also a small radial extra-acceleration of the form 

ApcTt = kHoVr- (1) 

In it A; is a positive numerical parameter of the order of unity to be deter- 
mined from the observations, Hq = (73.8 ± 2.4) km s~^ Mpc~^ = (7.47 it 
0.24) X 10~^^ yr~^ [44] is the Hubble parameter at the present epoch, defined 
in terms of the time- varying cosmological scaling factor S{t) as Hq = S/S , 

and Vr is the component of the velocity vector v of the test particle's proper 
motion about the central body along the common radial direction. The 
radial velocity for a Keplerian ellipse is [45] 



naesin/ 

where n is the Keplerian mean motion, a is the semi-major axis, and / is 
the true anomaly reckoning the instantaneous position of the test particle 
along its orbit: Vr vanishes for circular orbits. 
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The consequences of eq. ([T]) on the trajectory of the particle can be 
straightforwardly worked out with the standard Gauss equations for the 
variation of the Keplerian orbital elements [l5] which are valid for any kind 
of perturbing acceleration, whatever its physical origin may be. For the 
semi- major axis and the eccentricity they are 



da 
dt 



de 
dt 



{An,inf + AT [cos/ + i(l-^)]}- 



(3) 



In eq. ^ p = a(l — e^) is the semi-latus rectum, and Aji and At are the ra- 
dial and transverse components of the disturbing acceleration, respectively: 
m our case, eq. © is entirely radial. In a typical first-order perturbativB 
calculation like in the present case, the right-hand-sides of eq. ([3]) have to 
be computed onto the unperturbed Keplerian ellipse, characterized by 



P 



1 + e cos / 

and integrated over one orbital period by means of 

1 /r\2 1 



dt 



n 



fdf. 



(4) 



(5) 



It turns out that both the semi-major axis a and the eccentricity e of the 
test particle's orbit secularly increase according to 



(d) 



2kaHo (1 - Vl - e2 



(6) 



The formulas in eq. ([6]), which were obtained by taking an average over a 
full orbital revolution, are exact to all order in e. 

Since eMoon = 0.0647, it turns out that eq. ([6]) is able to reproduce the 
measured anomalous increase of the lunar orbit for 2.5 < A; < 5. Moreover, 
for such values of k eq. ([6]) yields an increase of the lunar semi-major axis 
of just 0.3 — 0.6 mm yr~^. It is, at present, undetectable, in agreement with 



'^Indeed, it can be easily inferred that eq. ([T]) is of the order of 10~^^ m s~^ for the 
Earth's motion around the Sun, while its Newtonian solar monopole term is as large 
as 10"^ m s"^ The same holds for the Earth-Moon system as well. Indeed, eq. ([T)) 
yields about 10^^^ m s^^ for the lunar geocentric orbit, while the Newtonian monopole 
acceleration due to the Earth is of the order of 10^^ m s^^. 
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the fact that, actually, no anomalous secular variations pertaining such an 
orbital element of the lunar orbit have been detected so far. If we assume the 
terrestrial semi-major axi^ a® = 1.5 x 10^^ cm as an approximate measure 
of the astronomical unit and consider that = 0.0167, eq. ([6]) and the 
previous values of k yield a secular increase of just a few cm yr~^. Also in 
this case, it can be concluded that eq. if applied to other situations for 
which accurate data exist, does not yield results in contrast with empirical 
determinations for a and e. Indeed, for the eccentricity of the Earth eq. 
§1, with 2.5 < A; < 5, yields (e) = (1.7 - 3.4) x lO"^^ yj.~i. Actually, 
such an anomalous effect cannot be detectable since, according to Table 3 
of Ref. |51] . the present-day formal, statistical accuracy in determining e 
from the observations amounts just to 3.6 x 10~^^; it is well known that the 
realistic uncertainty can be up to one order of magnitude larger. Similar 
considerations hold for the other planets. 



3 Conclusions 

Here we do not intend to speculate too much about possible viable physical 
mechanisms yielding the extra-acceleration o Heq. ©. 

It might be argued that, reasoning within a cosmological framework, 
the Hubble law may give eq. ([1]) for k = 1 ii the proper motion of the 
particle about the central mass is taken into account in addition to its purely 
cosmological recession which, instead, yields the well-known locaQ extra- 
acceleration of tidal type [46l [25l l47] 

^cosmoi = -qoHQr, (7) 

where go = — (^S/S^ Hq"^ is the deceleration parameter at the present 
epoch. 

On the other hand, our empirical results, which are not in contrast with 
other observational determinations for the Moon and the other planets of the 
solar system, may be simply interpreted, in a purely phenomenological way, 
in terms of a radial extra-acceleration proportional to the radial component 

^Actually, the astronomical unit is neither the semi-major axis of the Earth's orbit 
nor its average distance from the Sun (r) = a(l + e^/2) [3]. On the other hand, it is 
00 = 1.00000018 an ( |http:/7ssd.jpLnasa.gov/txt/p-elem_t2.txt[ ). 

''See Ref. [41] and his quasi-Newtonian dynamics, intermediate between Newtonian and 
modified Newtonian dynamics. 

*For a recent review of the infiuence of global cosmological expansion on the local 
dynamics and kinematics, see Ref. [12| . 
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Vr of the proper velocity of the test particle about its primary through a 
coefficient having the dimensions of T^^ and a magnitude close to that of 
Hq: its physical origin should not necessarily be of cosmological origin. 

Finally, we want to spend some words about the nature of the anomalies 
considered. As it was shown, the anomalous increase of the astronomical 
unit has attracted the attention of several researchers so far. It was deter- 
mined as a solve-for parameter by using different ephemerides which neither 
use the same dynamical force models nor the same observational records. 
Further processing of more extended data sets, with more accurate dynam- 
ical modeling, will be useful in shedding further light on such an anomaly. 
The authors of Ref. [l] conclude their review of the anomalous variation of 
the astronomical unit by writing at pag. 194: "If the reported increase holds 
up under further scrutiny and additional data analysis, it is indeed anoma- 
lous. Meanwhile it is prudent to remain skeptical of any real increase. In our 
opinion the anomalistic increase lies somewhere in the interval zero to 20 cm 
yr~^ , with a low probability that the reported increase is a statistical false 
alarm." On the other hand, it must be remarked that a clear definition for 
the change of the astronomical unit is still lacking since some researchers be- 
lieve that the astronomical unit is a redundant unit, like to the gravitational 
parameter GM of the Sun, which should, instead, be empirically determined 
from data processing as a solve-for parameter |48^ H9l [50] . It is likely that 
at the lAU meeting in 2012 a fixed numerical value for it will be adoptecH. 
Anyway, this would have the effect of just shifting the detected anomaly to 
another physical quantity. Concerning the Moon's orbit, the first report of 
the lunar eccentricity anomaly dates back to 2001 [6l; such a phenomenon 
is still here (Tj |8] , despite the increasing accuracy in LLR observations and 
modeling occurred in the last decade due to the steady efforts of a wide 
community of researchers engaged in LLR science and technology. On the 
other hand, it is not certainly unreasonable to expect that further modeling 
of classical effects occurring in the lunar interior may finally be able to ex- 
plain the observed anomaly. Anyway, until it will actually happen, looking 
for alternative explanations remains a task worth being pursued. 



References 

[1] Anderson J. D., Nieto M. M., 2010, Astrometric solar-system anoma- 
lies. In: Klioner S.A., Seidelmann P.K., Soffel M.H., (eds.) Relativity in 
Fundamental Astronomy: Dynamics, Reference Frames, and Data Anal- 

^W. M. Folkner, personal communication to the author, February 2011. 



5 



ysis, Proceedings lAU Symposium No. 261, Cambridge University Press, 
Cambridge, 2010, pp. 189-197 

[2] Krasinsky G. A., Brumberg V. A., 2004, Celest. Mech. and Dyn. Astron, 
90, 267 

[3] Standish E. M., 2005, The Astronomical Unit now. In: Km'tz D. W. 
(ed.). Transits of Venus: New Views of the Solar System and Galaxy, 
Proceedings of the lAU Colloquium No. 196, Cambridge University Press, 
Cambridge, 2005, pp. 163-179 

[4] Pitjeva E. V., Standish E. M., 2005, The Astronomical Unit now. In: 
Kurtz D. W. (ed.). Transits of Venus: New Views of the Solar System 
and Galaxy, Proceedings of the lAU Colloquium No. 196, Cambridge 
University Press, Cambridge, 2005, p. 177 

[5] E. V. Pitjeva, E. V., 2008, private communication to Noerdlinger P. 

[6] Williams J. G., Boggs D. H., Yoder C. F., Ratcliff J. T., Dickey J. O., 
2001, J. Geophys. Res., 106, 27933 

[7] Williams J. G., Dickey J. O., 2003, Lunar Geophysics, Geodesy, 
and Dynamics. In: Noomen R., Klosko S., Noll C, Pearl- 
man M. (eds.) Proceedings of the 13th International Workshop 
on Laser Ranging "Toward Millimeter Accuracy", Wshington 
D.C., October 7-11 2002. NASA/CP-2003-212248, p p. 75-86. 
http: / / cddis.nasa.gov/lwl3/ docs / papers / sci_williams_lm.pdf 

[8] Wilhams J. G., Boggs D. H., 2009, Lunar Core and Man- 
tle. What Does LLR See? In: Schilhak S. (eds.) Pro- 
ceedings of the 16th International Workshop on Laser Rang- 
ing "SLR-the next generation", Poznah, October 13-17 2008. 
[http:/ /cddis. gsfc.nasa.gov/lwl6/docs/papers/sci_l_Williams_p.pdf 

[9] Ni W.-T., 2005, Int. J. Mod. Phys. D, 14, 901 

[10] Kopeikin S. M., 2007, Relativistic Reference Frames for Astrometry and 
Navigation in the Solar System. In: Belbruno E. (ed.) NEW TRENDS 
IN ASTRODYNAMICS AND APPLICATIONS III. AIP Conference Pro- 
ceedings, Volume 886, pp. 268-283 

[11] Lammerzahl C, Dittus H., 2007, Int. J. Mod. Phys. D, 16, 2455 

[12] Carrera M., Giulini D., 2010, Rev. Mod. Phys., 82, 169 



6 



[13] Kopeikin S. M., 2010, Beyond the standard lAU framework. In: Klioner 
S. A., Seidelmann P. K., Soffel M. H. (eds.) Relativity in Fundamental 
Astronomy: Dynamics, Reference Frames, and Data Analysis, Proceed- 
ings of the International Astronomical Union, lAU Symposium, Volume 
261, Cambridge University Press, Cambridge, pp. 7-15 

[14] Li X., Chang Z., Li M., 2010. larXiv:1001.0066l 

[15] Li X., Chang Z., Mo X., 2010. larXiv:1001. 26671 

[16] Goldhaber A. S., Nieto M. M., 2010, Rev. Mod. Phys., 82, 939 

[17] Speliotopoulos A. D., 2010, Gen. Relativ. Gravit., 42, 1537 

[18] Zhang W. J., Li Z. B., Lei Y., 2010, Chinese Sci. Bull., 55, 4010 

[19] Sharma B. K., 2011, Earth, Moon, and Planets, 108, 15 

[20] Lammerzahl C, 2011, Testing Basic Laws of Gravitation- Are Our Pos- 
tulates on Dynamics and Gravitation Supported by Experimental Evi- 
dence? In: Blanchet L., Spallicci A., Whiting B. (eds.) Mass and Mo- 
tion in General Relativity Fundamental Theories of Physics, Volume 162, 
Springer, Berlin, pp. 25-65 

[21] Zhang W. J., Kelley N., 2011, Advanced Sci. Lett., 4, 574 

[22] lorio L., 2005, J. Cosmol. and Astropart. Phys., 09, 006 

[23] Mashhoon B., Singh D., 2006, Phys. Rev. D., 74, 124006 

[24] 0stvang D., 2007, Gravit. Cosmol, 13, 1 

[25] Mashhoon B., Mobed N., Singh D., 2007, Class. Quantum Gravit., 24, 
5031 

[26] Khokhlov D. L., 2007. larXiv:0710.5862l 
[27] Noerdlinger P., 2008. llarXiv:0801.3807l 

[28] Lammerzahl C, Preuss O., Dittus H., 2008, Is the Physics Within 
the Solar System Really Understood? In: Dittus H., Lammerzahl C, 
Turyshev S. G. (eds.) Lasers, Clocks and Drag-Free Control Exploration 
of Relativistic Gravity in Space, Astrophysics and Space Science Library 
Volume 349, Springer, Berlin, pp. 75-101 



7 



[29] Verbiest J. P. W., Bailes M., van Straten W., Hobbs G. B., Edwards 
R. T., Manchester R. N., Bhat N. D. R., Sarkissian J. M., Jacoby B. A., 
Kulkarni S. R., 2008, Astrophys. J., 679, 675 

[30] Lammerzahl C, 2008, Eur. Phys. J. Special Topics, 163, 255 

[31] Amin M. Y., 2009. larXiv:0912.2443l 

[32] Arakida H., 2009, New Astron., 14, 264 

[33] Miura T., Arakida H., Kasai M., Kuramata S., 2009, Publ. Astron. Soc. 
Jp., 61, 1247 

[34] Ito Y., 2009, Publ. Astron. Soc. Jp., 61, 1373 

[35] Li X., Chang Z., 2009. larXiv:0911. 18901 

[36] Lammerzahl C, 2009, Space Sci. Rev., 148, 501 

[37] Brumberg V.A., 2010, Celest. Mech. Dyn. Astron., 106, 209 

[38] Nyambuya G. G, 2010, Mon. Not. Roy. Astron. Soc, 403, 1381 

[39] Bel LI., 2010. larXiv:1005.5442l 

[40] Arakida H., 2010, Adv. Sp. Res., 45, 1007 

[41] Rasor N. S., 2010, Physics Essays, 23, 383 

[42] lorio L., 2011, Mon. Not. Roy. Astron. Soc, doi:10.1111/j.l365- 
2966.2011.18777.x 

[43] Arakida H., 2011, Gen. Relativ. Gravit., doi:10.1007/sl0714-011-1170- 
1 

[44] Riess A. G., Macri L., Casertano S., Lampeit H., Ferguson H. C., Fil- 
ippenko A. V., Jha S. W., Li W., Chornock R., Silverman J. M. , 2011, 
Astrophys. J., 730, 119 

[45] Brouwer D., Clemence G. M., 1961, Methods of celestial mechanics. 
Academic Press (New York) 

[46] Cooperstock F. I., Faraoni V., Vollick D. N., 1998, Astrophys. J., 503, 
61 



8 



[47] Klioner S. A., Soffel M. M., 2005, Refining the Relativistic Model for 
Gaia: Cosmological Effects in the BCRS In: Turon C, O'Flaherty K. S., 
Ferryman M. A. C. (eds.), Proceedings of the Symposium "The Three- 
Dimensional Universe with Gaia", October 4-7 2004, Paris. ESA SP-576, 
pp. 305-308 

[48] Klioner S. A., 2008, Astron. Astrophys., 478, 951 

[49] Fienga A., Manche H., Kuchynka P., Laskar J., Gastineau M., 2010, 
.arXiv: 1011. 4419] 

[50] Capitaine N., Guinot B., Klioner S., 2010, Proposal for the re-definition 
of the astronomical unit of length through a fixed relation to the SI me- 
tre. Rencontres de I'Observatoire Journees 2010 "Systemes de reference 
spatio-temporels" . New challenges for reference systems and numerical 
standards in astronomy. Observatoire de Paris Ecole Normale Superieure 
20-22 September 2010. p. 6 

[51] Pitjeva E. V., 2008, Use of optical and radio astrometric observations of 
planets, satellites and spacecraft for ephemeris astronomy. In: Jin W. J., 
Platais I., Ferryman M. A. C. (eds.) A Giant Step: from Milli- to Micro- 
arcsecond Astrometry Proceedings lAU Symposium No. 248. Cambridge 
University Press, Cambridge, pp. 20-22. 



9 



